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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 

RESEARCH  MEMORANDUM 

COMPARISON  OF  EXPERIMENTAL  WITH  CALCULATED  RESULTS  FOR 
THE  LIFTING  EFFECTIVENESS  OF  A FLEXIBIJE  45°  SWEPTBACK 
WING  OF  ASPECT  RATIO  6.0  AT  MACH  NUMBERS 
FROM  0.8  TO  1.3 
By  Richard  E.  Walters 

SUMMARY 


Tests  were  conducted  on  models  having  45°  sweptback  wings  with 
varying  degrees  of  flexibility  to  determine  the  effective  lift-curve 
slopes  in  order  to  evaluate  the  usefulness  of  a general  method  for  the 
prediction  of  the  effective  lift  ratio  by  a comparison  of  the  predicted 
values  with  the  test  results.  Tests  were  made  with  three  aspect-ratio-6.0 
taper-ratio-0.6  wings  having  45°  sweepback  of  the  q\oarter-chord  line  and 
NACA  65AOO9  streamwise  airfoil  sections.  The  Mach  number  range  covered 
was  from  0.8  to  1,5 • 


CxJTves  of  the  lift-curve  slope  and  pitch  damping  are  presented, 
effective  lift  ratio  ^Lo^^^Lar  determined  from  the  experimental 


The 


results  and  from  a representative  method  of  prediction  is  shown  as  a 

fmction  of  the  load  flexibility  parameter  Ct  qk  for  purposes  of  com- 

"ar 


parison.  A comparison  is  also  made  of  the  effect  of  different  assumed 
load  distributions  and  of  different  assumed  centers  of  pressure  on  the 
predicted  effective  lift  ratio. 


The  results  of  the  tests  showed  that  these  methods  predicted  values 
of  the  effective  lift-slope  ratio  which  were  within  5 percent  of  the 
experimental  values. 


INTRODUCTION 


With  the  increased  use  of  thin,  sweptback  wings  of  high  aspect  ratio, 
the  problem  of  elastic  deformation  has  assumed  primary  significance.  The 
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aerod^amic  characteristics  of  the  wing  can  no  longer  be  considered  inde- 
pendently of  the  structural  deflections  since  the  effect  of  wing  bending 
and  torsion  on  the  sectional  angle  of  attack  has  become  appreciable. 

There  have  been  many  and  varied  attempts  to  predict  and  evaluate 
the  change  in  the  lifting  effectiveness  of  a sweptback  wing  as  the  wing 
is  allowed  to  deflect  under  load  (for  example,  refs.  1 to  5)  • The  method 
of  reference  1 is  the  most  general  in  nature  and  allows  the  application 
of  both  arbitrary  load  distributions  and  wing  construction  to  the  problem. 
In  most  approaches,  however,  it  has  been  necessary  to  compromise  either 
the  structural  or  aerodynamic  aspects  in  order  to  obtain  a solution. 

(See  refs.  2 and  3.)  The  purpose  of  this  paper  is  to  compare  the  general 
methods  suggested  for  the  solution  with  the  experimental  results  of  flight 
tests  of  a representative  sweptback  wing  of  varying  degrees  of  flexibility. 

In  the  following  analysis  the  deformation  of  the  structure  is 
expressed  in  terms  of  a set  of  experimentally  determined  structural 
influence  coefficients.  Different  types  of  load  distributions  are  assimied 
and  applied  to  the  influence  coefficients  in  order  to  determine  the  pre- 
dicted effective  lift.  There  are  three  comparisons  to  be  made:  first, 

that  between  the  predicted  values  of  the  effective  lift  and  the  experi- 
mental resiilts;  second,  that  between  values  predicted  by  assuming  dif- 
ferent load  distributions  in  the  calculations;  and  thirdly,  that  between 
values  predicted  by  assuming  different  center-of-pressure  positions. 

The  experimental  results  were  determined  from  flight  tests  of  three 
rocket-powered  models  with  the  same  wing  plan  forms  but  varying  degrees 
of  wing  flexibility.  The  wings  were  of  aspect  ratio  6.0  and  taper 
ratio  0 .6,  and  had  NACA  65AOO9  free-stream  airfoil  sections . The  varia- 
tion in  wing  flexibility  was  due  to  the  differences  in  the  wing  inlays 
which  were  0.064-inch  Inconel,  0.032-inch  Inconel,  and  0.064-inch 
24s -T  aluminum  alloy. 

The  Mach  mmiber  range  covered  was  approximately  0.8  to  I.3  and  the 
Reynolds  mmiber  range  was  3*0  x 10°  to  8.0  x 10°  based  on  wing  mean  aero- 
dynamic chord.  The  models  were  flown  at  the  Langley  Pilotless  Aircraft 
Research  Station  at  Wallops  Island,  Va. 


SYMBOLS 


an  normal  acceleration,  g units 

b wing  span,  ft 

c chord,  ft 
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c 


mean  aerodynamic  chord 
lift  coefficient,  — 

^ qs 


lift-curve  slope,  57.5  per  radian 

da 

effective  lift-curve  slope  of  flexible  wing 
rigid-wing  lift-curve  slope 


section  lift -curve  slope 


■'m 


pitching-moment  coefficient. 


M 

qSc 


Cmq  + 


8 

k 


damping-in-pitch  coefficient,  + , per  radian 

bis.  bSS. 

2V  2V 

acceleration  due  to  gravity,  52.2  ft/sec'^ 


stiffness  parameter, 


e/L 


(e/L) 


ref 


L 

M 

N 

P 

q 

r 

S 


lift,  lb 

pitching  moment,  ft-lb 
normal  force,  lb 
load,  lb 

dynamic  pressxrre,  Ib/sq  ft 
body  radial  coordinate,  in. 
area  of  reference  panel,  sq.  ft 
wing  area,  sq  ft 
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torque,  in-lb 

weight  of  reference  panel,  lb 

total  configuration  weight,  lb 

longitudinal  body  coordinate,  in. 

spanwise  coordinate  normal  to  fuselage  center  line 

angle  of  attack,  deg 

local  angle  of  attack  of  flexible  wing,  deg 
local  angle  of  attack  of  rigid  wing,  deg 

change  in  angle  of  attack  caused  by  wing  deflection 

— — , radians  per  second 
3 dt 

angle  of  pitch,  radians;  angle  of  rotation  of  reference 
chord,  deg 

radians  per  second 


rotation  of  reference  chord  due  to  a unit  concentrated 
load  applied  at  reference  station 

distance  between  loading  axis  aind  assumed  center-of-pressure 
axis,  fraction  of  chord 

structural  influence  coefficients  for  angle -of -attack  change 
due  to  unit  concentrated  loads  applied  along  reference 
axis,  deg/lb 

structural  influence  coefficients  for  angle-of-attack  change 
due  to  unit  torque  applied  parallel  to  free  stream, 
deg/ ft-lb 

column  matrix 

square  matrix 
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H 

H 

LxJ 


diagonal  matrix 


diagonal  matrix  with  nonzero  elements  equal  to  1 


square  matrix  with  all  elements  equal  to  1 
row  matrix  with  all  elements  equal  to  1 


MODELS  AND  TESTS 
Models 


A sketch  of  the  models  tested  showing  the  pertinent  dimensions  is 
presented  in  figure  1.  Photographs  of  the  models  are  shown  in  figure  2. 
The  fuselage  was  a curved  body  of  revolution  with  a maximum  diameter  of 
6.77  inches  and  a fineness  ratio  of  10.  The  fuselage  ordinates  are  given 
in  table  I. 

The  three  models  tested  were  of  similar  construction  except  for  the 
metal  inlays  in  the  wings . The  wing  geometry  was  as  follows : aspect 

ratio  6.0,  taper  ratio  0.6,  free-stream  airfoil  section  NACA  65AOO9,  and 
45°  angle  of  sweep  of  the  quarter-chord  line.  The  wing  construction 
showing  the  inlays  is  presented  in  figvire  5(a).  The  different  inlays 
with  their  respective  9/l  values  were  as  follows;  0 .064-inch- thick 
Inconel  for  model  1,  with  G/L  of  -0.0075  degrees  per  pound;  0.052-inch- 
thick  Inconel  for  model  2,  with  G/L  of  -0.0112  degrees  per  pound;  and 
0 .064 -inch- thick  24  S-T  alianinum  alloy  for  model  5^  with  G/L  of 
-0.0224  degree  per  poxmd. 

Directional  stability  was  obtained  for  the  models  by  the  use  of  two 
vertical  tails  of  24  S-T  aluminm.  The  tail  plan  form  may  be  seen  in 
the  sketch  of  figure  1. 

The  models  were  equipped  with  four-channel  telemeters  which  provided 
measurements  of  normal  and  longitudinal  acceleration,  total  pressure,  and 
angle  of  attack. 

Dviring  the  coasting  portion  of  the  flight,  the  models  were  distiirbed 
in  pitch  by  successive  firing  of  eight  pvilse  rockets.  These  pulse  rockets 
were  located  in  the  fiaselage  in  groups  of  four  with  their  lines  of  thrust 
perpendicular  to  the  plane  of  the  wings.  (See  fig.  1.)  The  total  impulse 
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of  the  pvilse  rockets  used  was  about  6.1  pound-seconds  and  the  thrust- 
time cvirve  is  approximately  120  sin  39»3t  pound  from  zero  time  to 
0.08  second. 


TESTS 

Structural  influence  coefficients  were  measured  on  test  panels  which 
reproduced  the  wing  structure  as  closely  as  was  possible.  The  influence 
coefficients  for  the  angle-of-attack  change  were  determined  for 

loadings  on  the  25-percent-chord  line,  which  will  be  referred  to  here- 
after as  the  reference  axis  and  for  loadings  on  the  50 -percent-chord 
line.  A linear  variation  between  the  experimental  influence  coefficients 
obtained  along  the  25-  and  50 -percent-chord  lines  was  assumed  and  the 
torsional  coefficients  were  then  calculated  on  this  basis.  Figure  4 

shows  the  values  of  these  coefficients  for  model  3^  the  most  flexible 
of  the  series. 

The  models  were  launched  at  approximately  70°  from  the  horizontal 
by  means  of  a rail  launcher  (fig.  5)  • Model  propulsion  consisted  of  a 

65-lnch  HVAR  rocket  motor  as  a booster  with  a rocket  motor  sus- 

4 

tainer.  Atmospheric  data  were  determined  by  radiosonde  observations  and 
trajectory  and  flight  velocity  were  measured  by  an  SCR-584  radar  and  a 
CW  Doppler  radar  set,  respectively. 

The  variations  of  dynamic  pressure  and  Reynolds  number  with  Mach 
number  are  shown  in  figures  6 and  7* 


ANALYSIS 


Inasmuch  as  the  purpose  of  this  analysis  is  to  compare  the  results 
given  by  existing  methods  for  the  prediction  of  the  effect  of  wing  elas- 
ticity upon  the  rigid-wing  lift-curve  slope  with  experimental  values, 
it  would  be  well  to  state  the  primary  methods  which  have  been  suggested. 

There  are  two  general  approaches  to  the  problem  which  might  be 
termed  the  aerodynamic  and  the  structinral.  In  the  aerodynamic  approach, 
the  structural  aspects  are  not  \isually  developed  in  detail  and  their 
effects  are  accounted  for  through  assumed  deflection  curves  whereas  the 
main  emphasis  is  placed  on  the  aerodynamic  considerations.  Conversely, 
in  the  structural  approach,  the  aerodynamic  contribution  is  usually  esti- 
mated by  use  of  strip  theory  which,  in  some  cases,  includes  a so-called 
tip  correction,  and  the  structure  is  investigated  in  detail. 
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If  the  equation  representing  the  contribution  of  forces  and  moments 
at  one  station  to  the  change  in  angle  of  attack  at  another  station  is 
vrritten  as 


ttg  = 0pP  + ^ecP 


it  can  be  seen  that  the  structural  influence  is  manifested  directly 
through  the  influence  coefficients  0p  and  and  the  aerodynamic, 

through  the  loading  P and  the  moment  arm  ec.  Therefore,  the  tvo 
approaches  may  be  discussed  and  appreciated  through  an  investigation  of 
these  separate  quantities  and  and  P and  ec  along  with  their 

individual  effects  upon  the  lifting  effectiveness  of  the  wing. 


Aerodynamic  Approach 

For  this  investigation,  the  change  in  angle  of  attack  caused  by  the 
action  of  aerodynamic  and  inertia  loads  on  a flexible  wing  structure  is 
considered  to  be  the  sum  of  the  torsional  and  bending  contributions  of 
the  load  distribution.  The  twisting  effect  is  the  result  of  the  load 
distribution  having  a center-of-pressure  axis  displaced  from  the  refer- 
ence axis.  The  bending  of  the  wing  also  effectively  causes  a rotation 
of  the  free-stream  chord  because  the  wing  tends  to  bend  and  twist  normal 
to  the  reference  axis  which  is  swept  back  at  an  angle  to  the  free  stream. 

Various  methods  have  been  suggested  to  approximate  the  lift  distri- 
bution on  flexible  sweptback  wings.  Some  of  the  distributions  which  have 
either  been  used  in  the  previous  methods  or  appear  applicable  for  use 
are  as  follows ; Weissinger ' s simplified  lifting-surface  theory  (refs . 4 
and  5)  awd  empirical  methods  based  on  this  theory  or  on  lifting-line 
theory  (ref.  6)  for  subsonic  speeds;  linearized  lifting-svjrface  theory 
(ref.  2)  for  supersonic  speeds;  and  strip  theory  with  or  without  tip 
corrections  for  all  speeds,  as  used  in  reference  5 and.  the  calculations 
of  reference  1. 

The  effect  of  these  different  assumed  load  distributions  on  the 
effective  lift-curve  slope  ratio  which  is  the  ratio  of  the 

flexible-wing  lift-curve  slope  to  the  rigid-wing  lift-curve  slope,  can 
be  determined  by  evaluating  the  effective  lift  produced  by  these  distri- 
butions in  conjunction  with  experimentally  determined  influence  coeffi- 
cients. The  most  convenient  approach  to  the  problem  of  representing  the 
deflections  and  rotations  of  the  wing  structure  appears  to  be  in  the  form 
of  Influence  coefficients.  This  method  obviates  the  representation  of 
the  structural  deformation  as  a series  of  assumed  deflection  modes. 
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The  use  of  influence  coefficients  reduces  the  problem  to  the  solu- 
tion of  a set  of  simultaneous  equations;  this  procedure  is  facilitated 
by  matrix  notation. 


Structural  Approach 

Experimentally  determined  influence  coefficients  are  the  most  desir- 
able but  naturally  assume  the  wing  to  be  available  for  the  necessary 
testing.  When  this  is  not  the  case,  the  structural  behavior  of  the  wing 
must  be  approximated. 

The  most  frequently  applied  approximation  is  that  based  on  simple 
beam  theory  where  the  wing  is  assumed  to  be  cantilevered  at  an  "effective 
root."  The  effective  root  was  initially  considered  to  be  a line  normal 
to  the  elastic  axis  passing  through  the  intersection  of  the  elastic  axis 
and  the  fuselage  chord.  The  elastic  axis  is  usually  considered  to  be  the 
locus  of  the  section  shear  centers  without  consideration  of  the  effect  ■ 
of  root  restraint. 

This  method  does  not  adequately  represent  the  rigidity  of  the  tri- 
angular portion  of  the  wing  formed  by  the  wing  root  and  the  effective 
root.  A truer  representation  is  afforded  if  the  effective  root  is  moved 
outboard.  The  required  amount  of  movement  of  the  root  is  uncertain, 
inasmuch  as  the  exact  position  can  usually  be  determined  only  by  experi- 
mentation or  by  a very  detailed  analysis  of  the  wing  structure.  A good 
approximation  to  the  position  of  the  effective  root  is  that  formed  by  a 
line  normal  to  the  elastic  axis  passing  through  the  intersection  of  the 
fuselage  chord  and  the  wing  trailing  edge.  This  concept  is  more  fully 
explained  in  reference  'J } however,  the  use  of  influence  coefficients 
makes  the  consideration  of  the  elastic  axis  unnecessary. 

Approaches  assimiing  specific  deflection  curves  or  those  which  are 
based  on  geometric  or  struct\n:al  criteria  are  not  discussed  Inasmuch  as 
these  wings  are  so  constructed  as  to  fit  the  necessary  assumptions  and, 
consequently,  are  no  longer  "arbitrary"  wings. 


Method  Used  in  the  Present  Paper 

Assumptions .-  For  the  method  used  in  this  paper,  the  following 
assumptions  are  made: 

(l)  The  total  angle-of-attack  change  due  to  wing  flexibility  Ug 

is  a resixlt  of  wing  bending  and  torsion,  and  these  effects  may  be  sepa- 
rated and  treated  individually. 
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(2)  For  the  purpose  of  determining  the  load  distribution,  the  wing 
is  divided  into  a number  of  panels.  The  lift-curve  slope  of  the  panel 
is  assumed  to  be  that  of  the  chord  at  the  panel  midpoint.  The  load  on 
the  panel  is  assumed  to  be  a concentrated  load  acting  at  the  intersection 
of  the  panel  midchord  and  the  center-of-pressure  axis.  This  point  is 
called  a loading  point.  Figure  3(t>)  shows  the  division  of  the  wing  into 
the  reference-panel  areas  and  the  positions  of  the  loading  and  measuring 
stations  as  assmed  for  this  analysis. 

(3)  The  center-of-pressure  axis  is  assumed  to  be  at  a constant  per- 
cent of  the  chord.  This  assumption  is  maintained  throughout  the  investi- 
gation; however,  a means  of  treating  those  cases  in  which  the  center  of 
pressure  is  not  a constant  percent  of  the  chord  is  preser  ted  in  the 
"Analysis"  section. 

(4)  Aerodynamic  induction  effects  are  not  considered  after  the 
initial  load  distribution  has  been  assumed.  Strip  theory  is  used  to 
calculate  the  lift  caused  by  a structural  deformation  so  that  the  changes 
in  the  lift  on  a reference  station  do  not  influence  the  lift  on  any  other 
panel. 


Development  of  the  aeroelastic  equation.-  The  method  presented  here 
is  similar  to  that  of  reference  1 and  is  simplified  by  assuming  constant- 
chord  segments  and  not  using  integrating  matrices . The  basic  equation 
for  the  contribution  of  the  forces  and  moments  at  one  station  to  the 
change  in  angle  of  attack  at  the  same  or  another  station  is 


(1) 


where  P is  the  distributed  load  along  the  wing  span  and  T is  the 
twisting  moment  produced  by  the  displacement  of  P from  the  elastic  axis 
or  from  the  reference  axis  used  for  the  determination  of  0p. 


T = ecP 


therefore. 


Os  = + ec^)p 


(2) 


Since  the  structural  characteristics  of  the  wing  are  represented  as 
influence  coefficients,  the  loading  P must  be  expressed  in  a corre- 
sponding form.  If  the  Influence  coefficients  are  assumed  known  for  a 
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set  of  reference  stations,  the  effective  load  mtist  also  be  known  for 
these  stations  which  are  taken  at  the  odd  tenths  of  the  exposed  semispan. 
According  to  assumption  (2),  the  effective  load  on  a panel  can  be  assumed 
to  act  through  a point;  in  this  case,  the  point  is  the  reference  station. 
Simultaneous  considerations  of  all  the  loads  acting  on  all  the  stations 
resolves  the  problem  into  the  following  set  of  simultaneous  equations; 


^2.  ~ ^o(^P  ®*^^t)iO  * 

ttg^  = Pq(0p  + ■*■••• 

‘"35  = Po(i^P  + + * • * + "■  • • • 

Oq^  = Pq^9^  + + • • • 

“S9  = Pq(<^p  + ■*■••• 


> (5) 


The  subscripts  of  Ug  refer  to  the  spanwise  station,  in  tenths  of 

the  exposed  semispan,  at  which  the  change  in  angle  of  attack  is  measured. 
The  subscript  of  P refers  to  the  exposed  semispan  station,  also  in 
tenths  at  which  the  load  is  applied.  The  first  subscript  of  the  combined 
influence  coefficients  ^0p  + ec^rp^  refers  to  the  spanwise  station  at  which 

the  change  in  angle  of  attack  was  measured  and  the  second  that  at  which 
the  load  was  applied.  For  example,  + ec0Qi^^^  in  the  eqviation  for 


means  the  rotation  of  the  chord  is  measured  at  station  5 and  is  due 
^5 


to  a load  applied  at  station  5 


this  load  P,  to  a. 

5 ^=5 


The  amount  of  rotation  contributed  by 

+ ec0j)^5. 


,g^  is  equal  in  magnitude  to  P^(^P 


The  above  set  of  sim\iltaneous  equations  and  its  subsequent  manipu- 
lations may  be  most  readily  handled  by  matrix  notation. 

The  equations  leading  up  to  equation  (5)  when  rewritten  in  matrix 
form  become 


{S} ' * [® 


The  twisting  moment 


w 


ccmwEmm. 
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Hence, 


{“•s}  = [^J  + {p} 


(5) 


The  load  P on  any  reference  panel  equals  the  lift  on  the  panel 
minus  the  effect  of  wing  inertia 


Ppanel  - ^panel  “ ^npanel'^panel 


(6) 


where  is  in  g units,  and 


^ N Ltotal^ 
^ ^total 


(7) 


The  lift  on  an  elastically  deformed  wing  is 


ob/2 

L = 2q  / cc-  a dy 

J Q 


which  in  matrix  notation  is 


or 


L 


(8) 


In  equation  (8)  and  the  following  derivation,  the  matrix 


is  used  to  represent  the  different  types  of  loading  distributions  (sim- 
plified subsonic  lifting  surface  theory,  etc.)*  For  strip  theory 
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where  K represents  an  overall  reduction  In  lift.  If  this 


factor  is  ignored,  a procedure  which  is  justified  to  a certain  extent  at 
supersonic  speeds,  K Is  equal  to  one. 


If  eqimtions  (7)  and  (8)  are  substituted  into  equation  (6),  the 
following  expression  is  obtained: 


{p}  = 


r o n 
‘■a 


For  the  purpose  of  this  derivation,  the  last  term  in  this  equation 
is  in  an  inconvenient  form;  the  desired  form,  which  consists  of  a square 
matrix  multiplied  by  a col\jmn  matrix  of  the  local  angle-of-attack  values, 
can  be  obtained  by  a device  used  in  reference  8,  since 


which  becomes 


{p}  = <1=1^  P]  - 
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The  aeroelastic  equation  is  then  obtained  by  substituting  this  expression 
for  -[p}  into  equation  (5^ } the  resulting  equation  is 


lecj 


Structural  influence  coefficients  were  determined  experimentally. 
The  ratio  between  corresponding  influence  coefficients  of  the  reference 
wing  and  another  wing  of  similar  construction  and  plan  form  is  constant 
and  equal  to  the  reference  g/L  values  of  each  wing  where  the  parameter 
G/L  refers  to  a rotation  of  the  streaniwise  chord  at  a reference  station 
through  an  angle  G due  to  a concentrated  load  L applied  at  the  same 
or  a second  reference  station.  This  ratio  may  be  expressed  as 

— 9/1 

(0/L) reference 


where  (Q/i) reference  ^^^'ers  to  the  wing  for  which  the  influence  coef- 
ficients were  determined.  This  ratio  may  be  introduced  into  equation  (9) 
as  follows ; 


[w 


and  since,  by  definition, 


the  solution  of  the  equation  is 


Equation  (lO)  may  be  solved  in  a number  of  ways  with  Grout's  method  of 
reference  9;  probably  the  most  suitable  method  for  manual  computing 
machines . 
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If  the  center  of  pressure  is  located  at  a constant  percent  of  the 
chord,  the  matrix  jecj  in  equation  (10)  can  be  replaced  by 

ej^ij^  [_c],  and  if,  in  addition,  strip  theory  is  used  with  a factor  K 
equal  to  1,  equation  (lO)  becomes 

{a}  = [l]  - qCL^k  [cpp]  + [l]  - 2^  [zi] 

J (11) 

If  the  twisting  effect  is  neglected  this  equation  is  further  simplified 
to 

= [ij  - l0i^k|9p]  [i]  - a[i]  C3  fog].  (lla) 


and  if  the  inertia  effects  are  neglected  and  the  values  e are  assumed 
constant  along  the  span,  equation  (ll)  becomes 


[ce}  =.  [l]  - [cpp]  + Ci]  ^ (12) 


The  effective  lift-curve  slope  Ct_  niay  be  defined  as 


_ L 


■'a.  qSa„ 


where  ag  is  assumed  constant  along  the  span. 


Eqtxation  (8)  can  be  written  as 


I = Oi^qSog  = 2qC,  LlJ  K1  ^ {5} 
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so  that  the  effective  lift  ratio  is 


= 2 


(15) 


The  effects  of  substituting  a strip-theory  distribution  for  a more 
exact  distribution,  of  neglecting  the  effect  of  the  twisting  contribution, 
and  of  neglecting  the  effect  of  the  inertia  forces  are  investigated  in 
the  illustrative  example. 


ILLUSTRATIVE  EXAMPLE 


In  order  to  illustrate  the  differences  in  the  effective  lift  ratio 
Cl  |Ct  for  different  assumed  load  distributions  and  center-of-pressixre 
oe/  “r 

positions,  the  load-flexibility  parameter  qCx  k is  assumed  to  be  suf- 

ficiently  large  to  cause  a loss  of  about  forty  percent  in  the  lifting 
effectiveness  of  the  wing. 

The  experimental  influence  coefficients  for  the  wing  with  the 
0.06l4-inch  aluminum -alloy  inlays  are 


Measuring  Stations 


0 

1 

5 

5 

7 

9 


0 

0 

0 

0 

0 

0 

0 


Loading  Stations 


1 

5 

5 

7 

9 

0 

0 

0 

0 

0 

0 

-0.0006 

-0 .0011+ 

-0.0024 

-0 .0044 

0 

-0 .0008 

-0.0055 

-0 .0066 

-0.0124 

degrees 

0 

-0 .0015 

-0.001+5 

-0 .0096 

-0.0196 

pound 

0 

-0.0012 

-0 .001+6 

-0 .0106 

-0 .0222 

0 

-0 .0011 

-0.0041 

-0 .0104 

-0.0224 
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T 

0 

0 

0 

0 

0 

0 

0 

0.0005 

0.0012 

0.0021 

0 .0039 

0 

0 

0.0008 

0 .0035 

0.0065 

0 .0118 

degrees 

0 

0 

0 .0016 

0 .0047 

0 .0099 

0 .0203 

foot-pound 

0 

0 

0 .0014 

0 .0052 

0.0121 

0 .0252 

0 

0 

0.0014 

0.0051 

0 .0131 

0 .0281 

Equation.  (lO)  is  solved  by  using  a rigid-wing  lift  distribution  calculated 
from  reference  5*  Equation  (ll)  which  utilizes  strip  theory  is  also 
solved  and  the  different  results  are  compared  to  determine  the  effect  of 
the  load  distribution. 


The  expression  for 


number  of  0.8  is 


o -) 
‘•a 


Cl, 


according  to  reference  5 Q-  Mach 


r 0 -1 


0.842 

0 

0 

0 

0 

0 


0 

0,922 

0 

0 

0 

0 


0 

0 

1.06l 

0 

0 

0 


0 

0 

0 

1.155 

0 

0 


0 

0 

0 

0 

l.l6o 

0 


0 

0 

0 

0 

0 

0.920 
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For  this  wing  with  0.064-inch  aluminum-alloy  inlays,  the  weight  distri- 
bution is  given  by 


* 

* 

* 

0 .0345 

0 .0345 

0 .0345 

0.0301 

0.0301 

0.0301 

0.0260 

0.0260 

0.0260 

0 .0223 

0 .0223 

0 .0223 

0 .0188 

0 .0188 

0 .0188 

* 

* 

* 

0.0345 

0 .0345 

0 .0345 

0.0301 

0 .0301 

0 .0345 

0.0260 

0.0260 

0.0260 

0.0223 

0 .0223 

0 .0223 

0.0188 

0 .0188 

0 .0188 

(The  number  denoted  by  an  asterisk  is  luanecessary  since  the 
multiplying  elements  in  the  influence -coefficient  matrices  are  zero.) 
The  pertinent  geometric  characteristics  are  given  by 


0 = 


0.320 

0 

0 

0 

0 

0 

0.976 

0 

0 

0 

0 

0 


0 

0.381 

0 

0 

0 

0 

0 

0.912 

0 

0 

0 

0 


0 

0 

0.352 

0 

0 

0 

0 

0 

0.843 

0 

0 

0 


0 

0 

0 

0.323 

0 

0 

0 

0 

0 

0.774 

0 

0 


0 

0 

0 

0 

0.295 

0 

0 

0 

0 

0 

0.706 

0 


0 

0 

0 

0 

0 

0.270 

0 

0 

0 

0 

0 

0.638 


The  load-flexibility  parameter  CL^^^qk  is  assumed  to  be  10,000  and  the 
rigid-wing  angle  of  attack  Og  is  assumed  constant  at  1°  along  the  span. 
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The  center-of-pressure  axis  is  assumed  at  the  0.25-chord  line,  which 
is  the  reference  axis  0p. 

Solution  of  equation  (lO)  using  the  above  values  yields  the  result- 
ant load  distribution  Substituting  this  value  into  eqmtion  (13) 


yields  a value  of 


of  0.646  for  this  case. 


Changing  the  load  distribution  from  that  of  reference  5 to  a strip- 

theory  distribution,  for  the  same  conditions,  yields  a value  of  — ^ 

of  0.652,  a change  of  less  than  1 percent  for  this  case  of  heavy  wirig 
loading . 


If  the  second  example,  that  using  strip  theory,  is  changed  by 
assuming  the  center-of-pressure  axis  to  be  on  the  0.50-chord  line  rather 
than  the  0.25-chord  line,  a comparison  may  be  made  to  determine  the  effect 
of  center-of-pressure  position  for  this  case. 

Changing  only  the  value  of  e in  equation  (lO)  and  solving  the 

resulting  equation  gives  a value  of  -ise.  of  0.684  or  about  3 percent 

difference  from  the  previous  case  for  the  0.25-chord  loading. 


RESULTS  AND  DISCUSSION 


Figures  8 and  9 show  the  results  of  the  flight  tests  of  the  three 
models . The  procedure  used  in  reducing  the  data  from  the  telemeter 
records  and  radar  observations  to  the  final  forms  as  presented  in  the 
figuires  is  explained  in  reference  10. 


In  figure  8 the  experimental  values  of  Cj^  against  Mach  number 
are  presented.  Figure  10  shows  the  extrapolation  of  Cj;^  to  obtain  the 
rigid-wing  values.  Figirre  9 shows  the  variation  of  Cjjj^  + Cjjj^  with  Mach 

number,  all  models  having  the  same  center-of-gravity  positions.  The 
experimental  values  of  the  pitching-moment-curve  slope  are  not  presented 
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in  this  paper  because  the  method  of  determination  used  assumes  a linear 
slope  of  the  pitching-moment  angle-of-attack  curve  throughout  the  test 
range.  Previous  Investigations  have  shovm  the  pitching-moment  angle-of- 
attack  curve  for  this  wing  to  be  nonlinear  at  small  positive  and  negative 
values  of  the  angle  of  attack. 


The  curve  of  the  effective-lift  ratio 
flexibility  parameter  Cx^^q.k  is  shown  in  figure  11. 


against  the  load- 
Here  the  results 


of  a strip- theory  load  distribution  calculated  by  equations  (ll)  and  (13) 
of  the  analysis  section  have  been  presented  with  center-of-pressure  posi- 
tions of  0.25~percent  and  0.50-percent  chord.  If  these  two  center-of- 
pressure  positions  are  assmed  to  be  the  boundaries  of  the  forward  and 
rearward  center-of-pressure  movement,  then  most  of  the  experimental 
points  fall  within  these  two  limits.  Those  which  fall  outside  the  limits 
are  possibly  in  error  due  to  the  experimental  accuracy,  the  inability  to 
extrapolate  to  the  exact  rigid-wing  lift-curve  slopes,  or  a combination 
of  the  two  coupled  with  the  accepted  error  arising  from  the  approximate 
methods  used  to  calculate  the  limits. 


The  results  of  comparing  the  effective  lifting  characteristics  of 
this  sweptback  wing  as  predicted  by  approximate  methods  with  the  experi- 
mental values  should  not  be  accepted  for  all  sweptback  wings  unless 
allowances  are  made  for  the  aspect  ratio.  It  is  believed  that  the  results 
do  show  that  these  approximate  approaches  are  sufficient  to  predict  the 
flexible-wing  lift-curve  slope  of  wings  having  an  aspect  ratio  of  6 or 
higher.  Reference  11  shows  that  a similar  approach  which  uses  influence 
coefficients  and  strip  theory  predicts  the  elastic  lift  of  a 45°  swept- 
back wing  of  aspect  ratio  4 with  a good  degree  of  accuracy. 

The  differences  between  the  values  obtained  by  using  either  strip- 
theory  or  a more  exact  approximation  for  the  rigid-wing  lift  distribu- 
tion are  small,  about  1 percent  where  the  greatest  measured  loss  in  lift 
was  recorded.  Inasmuch  as  the  rigid-wing  lift-curve  slope  can  not  always 
be  determined  with  any  more  accuracy,  it  appears  that  strip  theory  would 
suffice  for  most  of  the  cases  encountered  in  practice. 

For  this  wing  of  45°  sweepback  and  aspect  ratio  6.0,  the  data  indi- 
cate that  the  primary  contribution  to  the  twist  of  the  wing  is  that  of 
bending.  When  the  center  of  pressure  is  assumed  to  be  the  50-percent- 
chord  line  rather  than  the  25-percent-chord  line  the  difference  is  only 
3 percent  for  the  worst  loading  condition.  In  practice,  it  is  unusual 
if  the  rigid-wing  lift-c\a:ve  slope  is  known  to  any  better  accuracy. 
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CONCLUDING  RiMAEKS 


The  resvilts  of  the  flight  tests  of  three  similar  models  with  varying 
degrees  of  wing  flexibility  have  been  presented.  These  experimental 
values  of  the  lift-curve  slopes  have  been  compared  with  the  results  cal- 
culated by  a general  method  for  the  prediction  of  the  lifting  effective- 
ness of  flexible  sweptback  wings.  The  effects  of  changes  in  the  assumed 
load  distributions  and  centers  of  pressure  have  been  investigated  also. 

The  results  of  the  comparisons  indicate  that  the  effect  of  changing 
either  the  load  distribution  or  the  center  of  pressure  is  small,  about 
1 percent  and  5 percent,  respectively,  for  the  cases  of  rigid-wing  lift 
distributions  investigated  in  the  illustrative  example. 

The  values  of  the  effective  lift-slope  ratio  as  predicted  by  an 
assumed  strip  theory  load  distribution  coupled  with  experimentally  deter- 
mined structural  influence  coefficients  show  good  agreement  with  the 
experimental  results.  The  agreement  between  the  predicted  and  experi- 
mental values  is  within  the  accuracy  with  which  the  rigid-wing  lift-curve 
slope  can  usually  be  determined  in  practice. 
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TABLE  I 

BODY  COORDINAIES 


X,  In. 

r,  in. 

0 

0 

1.00 

.342 

2.00 

.578 

4.00 

.964 

6.00 

1.290 

8.00 

1.577 

12.00 

2.074 

16.00 

2.472 

20.00 

2.775 

22 .00 

2.892 

22.75 

2.953 

24.00 

2.995 

28.00 

3.146 

32.00 

5.250 

56.00 

5.314 

J+0.00 

5.354 

44 .00 

3.504 

48.00 

5.219 

52.00 

3.074 

56.00 

2.815 

60.00 

2.658 

64.00 

2.450 

66.70 

2.305 

67.70 

2.250 
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Figure  2.-  Concluded 


Sketch  of  model  wing  showing  wing  construction, 
construction  and  locations  of  loading  and  measuring  stations 
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S-.0200 

V. 

g> 

'^-.0/30 


Loading  station 

o 0.9^ 

Z nKpoitd 


° .7 
o .5 
.3 


0 .2  .4  .6  <5  J.0 

^posed  ^<t mi  span 

(a)  Influence  coefficients  due  to  bending  of  reference  axis. 


g^tKposdd  s<zmispin 

(b)  Influence  coefficients  due  to  twist  about  the  reference  axis. 

Figure  1+.-  Experimentally  determined  structural  Influence  coefficients 

for  model  3» 
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Variation  of  Ct  with  Mach  number  for  models  tested 
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Figiire  11.-  Comparison  of  experimental  resxilts  of  effective  lift-slope 
ratio  with  those  predicted  by  equations  (U)  and  (l5)» 


